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ABSTRACT 

A mathematical model of a vehicle being brought to rest in 
such a manner as to cause a skid is devised and simulated on an 
analog computer. This model includes the effects of the brake 
line pressure upon the braking of the vehicle and also the effects 
of the coefficients of friction between the tire and road surface 
and between the brake shoe and brake drum. Several controls to 
correct for the occurrence of a skid are simulated on the analog 
computer and a detailed analysis of the braking system under the 
influence of each control is performed. Lastly, a scheme for 
implementing the control which minimizes the stop time and the 
degree of skid is presented. 
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1. INTRODUCTION 



Skidding of a vehicle due to an excessive application of 

brake force for the road conditions prevailing causes many 

thousands of accidents every year* Rapid and heavy application 

of brakes in a panic stop situation often causes wheel lock up 

and subsequent loss of control of the vehicle. This also occurs 

in wet or icy road conditions. 

The minimum stopping distance possible for a vehicle is 

determined by the coefficient of friction between the tire and 

the road surface. When a vehicle is just on the verge of a skid 

the frictional force between the tire and road surface is at its 

greatest value, and consequently the deceleration of the vehicle 

is at a maximum, i.e., equal to [1> g, where is the coefficient 

s s 

of static friction between the tire and the road surface and g 

is the deceleration of gravity (see Appendix A for a proof of 

K 

this statement) . When skidding occurs the tire slides along the 
road surface. This causes a reduction in the coefficient of 
friction and a corresponding reduction in the frictional force 
between the tire and road surface, thus decreasing the magnitude 
of the vehicular deceleration. The optimum braking situation 
would thus be attained if the vehicle could be kept just at the 
verge of a skid. If a skid could be detected exactly at the 
instant of occurrence and a compensating force were automatically 
applied to eliminate the skid before wheel lock up occurred, then 
nearly ideal braking conditions could be realized under all road 
conditions . 
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Attempts have been made previously to develop "anti-skid" 
devices, and some such devices are presently in use on certain 
aircraft and many trucks. Previous efforts in this direction 
include the Dunlop System of England, the Maxaret Unit, the Anti 
Skid Device, the Westinghouse Decelostat Controller, the Hydro- 
Aire Hytrol system, and the Anti-Skid Braking system of Russia's 
Tu - 104 Jet Liner. [l] These devices all work satisfactorily to 
a greater or lesser extent, but most have the drawback that 
operation depends upon the sensing of either a complete wheel 
lock up or the attainment of a pre-determined rate of wheel 
deceleration, at which time a bang-bang type of control is initi- 
ated. Additionally, many of these devices require expensive, 
heavy installations, rendering them largely unsuitable for eco- 
nomical installation and operation on commercial vehicles. 

This thesis presents an analysis and evaluation of several 

electrically actuated alternate proposals designed to automatically 

compensate for an overbraking situation in order to minimize the 

degree of skidding and more nearly optimize the braking of a vehicle. 

A mathematical model for the system was developed and simulated on 

an analog computer, and three types of skid correcting controls 

were devised and tested. One form of control which linearly 

decreases the net force applied to the brake shoe until the skid 

» 

is eliminated and then permits a linear increase of the net braking 
force until a skid once again occurs was tested. This is called 
the double ramp control. Two forms of bang-bang controls were also 
tested. Finally a triple ramp control which continually "samples" 
the road surface to determine the optimum brake force was devised 
and tested. 
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Complete data were obtained for each type of control and a 
detailed analysis was made of the behavior of the braking system 
under the influence of each control. 

Lastly, a scheme for implementing the triple ramp control 
(which proved to be superior to any others devised) is advanced, 
although the physical system was not actually built. 
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2. DESCRIPTION OF THE PROBLEM 



For the purposes of this paper, a skid will be defined as 
that condition which exists when the actual linear velocity of 
the frame of the vehicle in question is greater than the product 
of the angular velocity of the wheel and the radius of the wheel* 
In terms of the symbology used throughout this paper 





e = V 


- RjCO, 


(i) 




e 0 


when v^>R^Cc) 


(2) 


and 


CD 

1! 

O 


when v = R^tO , 


(3) 



where v represents the linear velocity of the frame, R^ represents 
the radius of the wheel and oj is the angular velocity of the wheel* 

The term e corresponds to an error signal, the magnitude of which 
is a measure of the degree of skidding present. The product R^cO 
will be referred to as the velocity measured by the rotation of the 
wheel* When e = 0 there is no skidding and when e 0 the vehicle 
is moving faster than a point on the rim of the wheel, i.e„, the 
wheel is tending to lock up and a skid is occurring. The condition 
when RjCj>V would occur when the wheels are spinning, as, for example, 
when trying to accelerate on an icy road. 

A proposed means of sensing the linear velocity of the vehicular 
frame is to integrate the output from a linear accelerometer attached 
to the axle housing of the vehicle (to minimize effects of vehicle 
tilt during braking) , and to compare the integrator output with the 
velocity measured by the rotation of the wheel. The initial condition 
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input to the integrator would be the velocity measured by the 
rotation of the wheel at the instant of application of the brake 



force 



When an error signal occurs, indicating skidding is present, 
an automatic corrective force (F ) would be generated to reduce 
the total force F applied to the brakes. represents the force 

applied by the vehicle operator to the brakes and is considered 
constant for all subsequent analyses. 

These concepts are illustrated in Fig. 1 where a block diagram 
of the proposed system is shown^ 6 ^ 



Fig. 1. Block diagram of the proposed braking system. 

From elementary physics it is known that the coefficient of 
friction between two surfaces is dependent upon the relative velocity 
between the surfaces. A sketch of the behavior of the coefficient 
of friction as a function of the relative velocity between the two 
surfaces is shown in Fig. 2. 



Vo “■ R i 






Y 



Fc 
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Maxinnum value of coefficient 
of f irict ion 





RELATIVE ^ 

VELOC IT Y 



Fig. 2. Coefficient of friction as a function of relative velocity. 



Just before relative motion begins, the coefficient of friction 
between two surfaces is at its maximum value. Prior to reaching 
this maximum value, the coefficient of friction is a variable that 
adjusts itself in such a manner so that the force tending to cause 
a relative motion is just compensated by the reactive force of 
friction (this is before any relative motion takes place). After 
sliding motion starts (as between a brake shoe and drum, or between 
a tire and road surface for a skidding situation) the coefficient of 
friction varies as shown. 

A piecewise linear approximation to this behavior is shown in 
Fig. 3. 




Fig. 3. Piecewise linear approximation to the behavior 
of the coefficient of friction. 
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The term corresponds to the coefficient of static friction 
which is a maximum just before relative motion begins, and which is 
always greater than the coefficient of kinetic or sliding friction. 
The equation for this piecewise linear approximation is 



jUL— } for O < V V c 



(4) 



yU= — k/ (v- Vc) ; -for v c ^ V4(Vc +• £ 2 . 



(5) 



The above relationship will be used to represent the coefficient 

of friction between the brake drum and brake shoe of the moving vehicle. 

The break point, v^, for the coefficient of friction between the tire 

and the road surface will be assumed to be high enough so that for 

normal speeds the coefficient of kinetic friction is some constant 

\i . As will be seen, this is not a restrictive assumption, since the 
K 

purpose of the corrective forces to be applied is to minimize the degree 
of skid experienced. Thus the relative velocity between the tire and 
road surface would be small enough to warrant considering the coefficient 
of kinetic friction to be a constant. 

Consider the idealized brake drum and shoe shown in Fig. 4. 




BRAKE 

shoe: 

F= Fin — Fc 



T — A/et Torque 

£,■= Radius of Wheel (Tike) 
Re~ Radius of BraKe Drum 



Ff 



Fig. 4. Diagram of brake drum and shoe. 
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In Fig. 4 the wheel is assumed to be rotating in a clock- 
wise direction with angular velocity u> and angular acceleration 
(D with the positive sense as shown. F represents the total (net) 
force acting on the brake shoe, and is the difference between F. 
(applied by the brake pedal) and F c (the corrective force). F^ 
is the force of friction, the value of which depends upon whether 
a skid is present. The equations describing this system are pre- 
sented below. 

( 6 ) 

(7) 

( 8 ) 

> f or e>0 (9) 

where m = the weight of the vehicle acting through the axis of 

2 

the wheel divided by the acceleration of gravity g (= 16 ft. /sec. ), 



F = F. - F 
in c 

Net torque = T = Fp. - R^F^ 
F f = f — m V ; for e = o 



h 



_ f /U-o- K' Ra (co ~ CO c) j-for co > C0 c 



( 10 ) 

/U 0 j for co ^ 6c ) c U 1 ) 

i.e., |Jb is the coefficient of friction between the brake drum and 



brake shoe (Fig. 3). Also, |i is the coefficient of friction between 
the tire and the road when a skid is present. 

Let K | R ^ — K | , so that Eq. (10) becomes 






__ r ^/u 0 - k, (co-coc) , J r or co^ co c 
( /OLo ) -for co cOo 

and Eq. (7) becomes 

T = R e F/<. -R, F f = - Ju 

where J = the moment of inertia of the wheel. 



( 12 ) 

(13) 

(14) 
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Consider the case where e = 0, i.e., skidding is absent, so 



Rf = - m v i j ( 15 > 

with the equality holding when the system is just at the verge of 
a skido The equation of motion, (14), becomes 

— J” CO — R , hn V + R £ F ft- . (16) 

Since e = 0, then 

V - R, CO . (17) 

Substituting (17) into (16) yields 

- Jo> = Rf m £ + Re F/i . (is) 

Solving (18) for <J0 yields 

&> = - R* F/4 . <19) 

J+ Rf'm 

When e > 0, i P e., skidding is present, and 

FFf “ ) where ; (20) 

then (14) becomes 

- J, co ~ Re F/t - R , K rn gf ; (21) 



SO 



a; = R./“* 171 $ — Re. . 
J" J 



(22) 



Note that for each of the relations for (JO the value of p. is 



a function of w. In the case where (JO < co then (JL = p. , and the 

c o 

equations above apply with p> replaced by p.^. When uo ^ O)^, then 



yCC K, Cco ~ CJc") * (23) 

The complete set of equations is indicated below: 

When e = 0, CO CO c l 

CO = ~_Rs_F Ao • (24) 

X+J^rn 
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When e > 0, CO <C. CO c i 

CO — R, Ak mq 

7 ~ + 



BaE/Ao 

J 



(25) 



When e = 0, CO ^ COc l 



cb - R £ Fk 



T+Rfm 



J_ co 



Re FK a . 
J +- R , a m 



(26) 



When e > 0 , CO ^ COc • 



CO - R a F K , ca - Ra P Ke 4- R, /x-k I t\ (j # (27) 

j 7 j 

In the last two relations, 

Ka = /J-o + K , COc . 

The analog computer simulation of this system of equations is 
shown in Fig. 5. Note that Appendix B contains a brief description 
of the analog computer symbols used in this paper and an explanation 

of the logic elements used to control various portions of the simu- 

. C2,3,9,5j 

lation. 

To determine the equations governing the forward motion of the 
vehicle, recall that the maximum possible deceleration for a given 
coefficient of friction p, ? between the tire and road surface is 



_ X hrv<*x — /•* 3 ) 

where fx - JUL s or ^/U ^ • 

Thus, when e = 0, i.e., no skidding, 

- x = - R , co ^ /<- s 3 j 



(28) 



(29) 



where p> = the coefficient of static friction between the tire and 
s 



the road surface. Also, before a skid occurs, 
X =r R , CO - v . 



(30) 
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Fig. 5# Analog computer simulation of brake drum and shoe. 



After a skid occurs, 

X = -/U K g ; (31) 

where = the coefficient of kinetic friction between the tire 
K 

and road surface which is considered constant for this simulation, 
and /X k < yU s . 

Thus, when e > 0, 

v = k g t + Vo ; 

where is the velocity of the vehicle when the skid begins. 

The analog computer implementation of these equations and 
ideas is presented in Fig. 6, which provides a simulation for V 
under conditions of skidding or skidless braking. 




Fig. 6. Analog computer simulation of vehicular velocity. 



The composite analog computer set up for both the brake drum 
and shoe and the vehicular velocity is shown in Fig. 7. This is 
intended to represent a conventional braking system in existence 
on present day vehicles and no corrective force is included in 
this diagram. 
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3. PROPOSED SKID CORRECTING CONTROLS 



Any form of corrective control must be initiated by the 
occurrence of a skid, which in turn is indicated by the presence 
of a positive error signal. The purpose of the corrective control 
must be to lessen the net force F applied to the brake shoe when 
a skid is sensed. When the skid has ceased, the corrective force 
F^ must decrease to permit a reapplication of the full braking 
force to the brake shoe. Thus the value of F c will depend upon 
whether a skid is present or not. 

The first form of control to be tried is called a double ramp 
control, and a typical shape for the control is shown in Fig. 8. 




An examination of Fig. 8 indicates that the maximum value 

attained by F^ is F^ (the force applied by the vehicle operator) . 

Since F = F . - F , and since there can never physically be a 

me v j j 

negative force on the brake shoe, F - F. . The value of K. is 

’em 3 

not necessarily equal to K^. 

An analog computer simulation of the double ramp control is 
shown in Fig. 9. 
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Fig. 9. Analog computer simulation of the double ramp control. 



The integrator indicated in Fig. 9 may be in the compute mode 
for the entire problem simulation, due to the presence of the feed- 
back limiter at the output of the integrator which prevents F^ 
from ever going negative. 

A bang-bang control was tried for the second form of corrective 
force. Since any physically realizable bang-bang control would have 
either a dead zone or a time delay, one control with a dead zone 
and one with a time delay were tried. The expression "dead zone" 
is meant to imply that there is some positive value of error volt- 
age above which the corrective force is full on until the error is 
once again reduced to zero. The analog computer implementation of 
the bang-bang control with a dead zone is shown in Fig. 10, where Oi 
is the value of the dead zone voltage. 




comparator 





OI 




F c 



Fig. 10 o Analog computer implementation of the 
bang-bang control with a dead zone. 
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The expression ’’time delay” as used here implies that the 
error exists for some time before corrective action is initiated. 

In order to simulate the time delay, a Wood fourth-order approximant [2] 
time lag circuit was set up to delay the error voltage by one second 
(real time delay), which in problem time corresponded to a delay time 
of 0.1 second. The analog computer diagram for the bang-bang control 
with a real time delay, with the potentiometer settings needed to 
produce a one second real time delay, is shown in Fig. 11. 




Fin 



Fig. 11. 



e(t-i)>o 

r" 

> 

Fc 

Analog computer implementation of the 
bang-bang control with a real time delay. 



The third form of control to be simulated is called the triple 
ramp control, and a typical plot of corrective force F^ versus time 
plus the complementary plot of F (= F^ - F^) versus time is shown 
in Fig. 12. 
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Fig. 12. Triple ramp control. 

The triple ramp control is of the following form. F^ = 0 
until a skid occurs. When a skid occurs, e > 0 and F^ begins 
rising rapidly with a slope to a maximum value of F^ or until 
e = 0. When e = 0, F^ then begins decreasing with a slope equal 
to - until e > 0 once more. Up to this point the control is 

identical to the double ramp control. However, the second time 
that e > 0 a separate monitoring device with a "memory" (a Mylar 
capacitor, for example) holds the voltage at the value corresponding 
to the corrective force at which the skid occurred. 

At the instant when e > 0 for the second time in Fig. 12, F^ 
would once again increase rapidly with slope until e = 0. Then 
F would again begin decreasing with slope - K. until it is just 

C 4 

above the corrective force at which skidding occurred previously. 
Then F^ would be allowed to decrease more slowly with a slope - 
(thus permitting a gradual increase in F until a skid once again 
occurs). The cycle then repeats itself. 
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The design of this particular type of control was motivated 
by the considerations expressed in the Introduction, i.e., that 
the optimum braking situation would be realized if the vehicle could 
be kept just at the verge of a skid. The storage of the value of 
F at which a skid occurs endows the control with a form of "intel- 



ligence" which permits the net force applied to the brake drum to 

be just less than the optimum force. The gradual increase in F 

permitted by the third portion of the control (with a slope of - K ) 

is necessary to permit the control to adapt to improving (i.e., less 

slippery) road conditions. The first cycle of the control, where 

F^ increases rapidly with a slope of and then decreases with a 

slope of - until the skid occurs once again, is necessary to 

permit the control to "test" the road surface for the net force 

required to just cause a skid. The analog computer implementation 

of this control is shown in Fig. 13. 

Integrator No. 1 provides the ramp with a slope of K^, - K^, 

or - as appropriate. This integrator is in the compute mode for 

the entire simulation time and F^ is kept from going negative by the 

lower bound of the feedback limiter. The upper bound on the limiter, 

and consequently on F^ , is determined by F . Potentiometer No. 1, 

labeled C.R./10, determines the value of F at which the slope 

c 

changes from - to - ; the ratio of F^^/F^^ (see Fig. 12) 

is termed the Correction Ratio (C.R.). It is this ratio which 

determines the setting of potentiometer No. 1. When integrator No. 2 

is in the track (reset) mode switch No. 1 is closed (FF = 1) and we 

desire the slope of F to be either K. or - K, (but not - K ), since 

c 3 4 c 

this is the time period between the first and second occurrences of 
the positive error signal. Thus to positively prevent relay No. 1 
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Fig. 13. Analog computer implementation of the triple 
ramp control. 
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from throwing during this phase of the cycle potentiometer No. 2 
is set lower than potentiometer No. 3 to insure that C = 1 and 
C = 0, where 1 and 0 correspond to a logical true and false, 
respectively . 

The flip flop (the output of which is denoted by FF and FF) 
is an R - S flip flop, and is used to control the mode of integrator 
No. 2 as well as the state of switches No. 1 and 2. 
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4 . THE UNCONTROLLED SITUATION 



The analog computer simulation of the uncorrected braking 
of the vehicle (shown in Fig. 7) was set up on the COMCOR Ci 
5000 Analog Computer, using the following arbitrary parameters: 



Thus 

and 

where 



K 1 = 0.01 

uo =60 rad. /sec. 
c 

|i =0.8 

o 

= 0.8 - 0.01 (uo - 60), (JO ^ U) 

\i = 0.8, uo < u) 

\J* = coefficient of friction between the brake drum 
and shoe. 

a = 0.75 => static friction between tire and road surface, 
s 



H K - 0.5 



kinetic friction between tire and road surface. 



V =75 ft. /sec. = initial vehicular velocity, 
o 



uo =75 rad/sec. 
o 



initial rotational velocity of the wheel. 



R 1 = 1 ft. 

R 2 = .6 ft. 

2 

J =25 ft. lb. sec. 

M = 93.75 slugs 

(Weight = mg = 2000 lbs. = weight of vehicle acting 
through axis of wheel). 

A 7000 lb. input force (corresponding to F^) was simulated as 
a 70 volt step input occurring at time = 0 sec. 

Six different variables were plotted for each run on a Brush 
Instruments Mark 200 (6 channel) Recorder. They were 

cu —> the angular acceleration of the wheel, in rad/sec, 

(JO => the angular velocity of the wheel, in rad/sec. 



29 



the actual linear velocity of the vehicle, 
in ft/sec. 

F the net force applied to the brake shoe 

(= F. - F ) , in lbs . 
m c 

e — > the error signal (= v - R^uo) , in ft/sec. 

F the corrective force, in lbs. 

c 

The uncorrected system response of the vehicle, i.e., the 
response obtained when no corrective force is applied, is shown in 
Fig. 14. From examination of the uncorrected system plot it is 
seen that the wheel locks in 0.96 sec. and skidding occurs 0.16 
sec. after application of the braking force. The skid is present 
until the vehicle comes to rest 4.64 sec. after the force is initially 
applied. 

Measuring the slope of the velocity output plot gives an 

2 

average deceleration of 16 ft/sec. (excluding the first 0.1 sec. 

of the problem, where skidding does not occur). This is expected, 

2 

since with p = 0.5 the deceleration during a skid is 16 ft/sec. 

K. 

(= 0.5 x 32 ft/sec. ^). 

To completely understand the shapes of the various plots 
displayed in Fig. 14, recall that the coefficient of friction between 
the brake drum and brake shoe is of the general shape indicated in 
Fig. 3, with a break point at = 60 rad/sec. The initial rotational 
velocity of the wheel was od^ = 75 rad/sec., and thus at the time of 
simulated application of the brakes the coefficient of friction |i> 
between the brake drum is 

jj. = 0.8 - 0.01 (75 - 60) = 0.65. 

The value of p will increase as the relative velocity between 
the brake shoe and brake drum decreases, until it has attained a 
maximum value of 0.8. After the problem begins there is a time lapse 
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Fig, 14, Unoorr acted braking situation. 
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of 0.16 sec. before skidding starts. During this initial 0.16 
sec. p> is increasing since the vehicle is slowing down* However, 
p> is sufficiently small during this time interval so that the net 
torque T (= “ ^l^P acting on the wheel is not sufficient to 

cause a skid. When the velocity has decreased to 70 miles/hr., 
and (JO is down to 70 rad/sec. (since a skid has not occurred yet), 
the value for p, has increased to 

p. = 0.8 - 0.01 (70 - 60) = 0.7. 

This value for p> is sufficiently large to cause the net torque 
to be just equal to that value necessary to cause a skid. At the 
instant a skid occurs the coefficient of friction between the tire 
and road surface drops abruptly from p> = 0.75 to p = 0.5. F also 

S K t 

decreases abruptly due to the occurrence of the skid and the net 

torque increases rapidly, causing the wheel to decelerate rapidly 

and lock after 0.96 sec. total elapsed time. 

Note that the plot for the radial acceleration of the wheel UO 

indicates that at time = 0 + the acceleration drops sharply to a 

2 

negative value (- 23 rad/sec. ) and then slowly becomes more negative 
until the time when a skid occurs. During this time the wheel is 
still slowing down gradually, and the deceleration (or negative 
acceleration) is increasing slowly as the value of p (between the 

brake shoe and brake drum) increases. At the point when the value 

• 2 

of UO becomes just less than - 24 rad/sec. , the skid occurs, and (JO 

2 

abruptly becomes much more negative (- 76 rad/sec. ) and decreases 

2 

rapidly until it reaches a constant - 94 rad/sec. The significant 

factor to be noted here is that the skid occurs just when W = 

2 ‘ 2 
- 24 rad/sec. . This corresponds to a value of V = - 24 ft/sec, 

(since = 1 ft., and a skid has not occurred yet). This is in 
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agreement with the previous statement that the maximum deceleration 

possible without a skid occurring is 

2 2 

8 = - 0.75 x 32 ft/sec v = - 24 ft/sec. . 
s 

• 2 

The leveling off of the ^ plot at - 94 rad/sec. occurs at 
the time when the rotational velocity w crosses the 60 rad/sec. 
mark which is the value established for the critical value of 
rotational velocity below which M* becomes a constant = 0.8. 

Since |l is constant and 

F f = ^mg 

then the net torque T is 

T = - R^p^mg = a constant, and (JO is thus a 

constant „ 

The magnitude of the error rises to a maximum value and then 
decreases to zero when the vehicle has completely stopped. This 
behavior is expected when it is noted that the error is merely the 
difference between the linear velocity and the rotational velocity 
of the wheel. 

One other comment concerning the response of the uncorrected 
system is appropriate. Inspection of the vehicular velocity v plot 
reveals that a steeper slope is present on the plot prior to the 
beginning of a skid. After the skid begins the slope becomes less 
steep and remains constant until the vehicle comes to rest. The 
deceleration of the vehicle is thus greater before the skid begins 
than it is after skidding occurs, which is in accordance with the 
expected results. 

The important fact to notice is that during a large part of 
the path the vehicle is in what would be an uncontrolled skid, as 
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shown by the large magnitude of the error. This uncontrolled skid 
is a matter of as much concern as the lengthened time to stop the 
vehicle. This lengthened time is caused by the effect of the 
smaller coefficient of friction which acts between the tire and 
the road surface when skidding is occurring. The optimum time of 
stopping, assuming braking takes place just at the verge of skidding, 
is 

75 ft/sec. ^ ^ 

t optimum = — - a =3.13 sec., 

2 

where 24 ft/sec. is the maximum possible deceleration using the 
values assumed for this simulation. 
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5. THE DOUBLE RAMP CONTROL 



The various proposed controls described in Section 3 were 

implemented in an attempt to obtain a more nearly optimal braking 

situation. The first of these corrective forces was the double 

ramp control shown in Fig. 8. A total of 21 separate runs were 

made with this form of control. was initially held constant 

at 40000 lbs/sec., and was varied from 5000 lbs/sec. up to 

60000 lbs/sec. Then was held constant at 60000 lbs/sec. and 

varied from 5000 lbs/sec. to 50000 lbs/sec. Finally, two 

runs were made with K 0 and K. set at 80000 lbs/sec. and at 100000 

3 4 

lbs/sec. Representative samples of the behavior to be expected 
from this form of control are shown in Fig. 15 where = 60000 
lbs/sec. and = 35000 lbs/sec., and in Fig. 16 where = 

100000 lbs/sec. and = 100000 lbs/sec. 

An examination of Fig. 15 indicates that for the specific values 
of and chosen (i.e., 60000 lbs/sec. and 35000 lbs/sec. respec- 
tively) the stop time has been increased to 5.47 sec., whereas the 
stop time under controlled conditions was only 4.64 sec. However, 
the total time during which a skid is occurring has been reduced to 
2.08 sec. and the magnitude of the skid has been reduced to approx- 
imately 2 ft/sec. during the entire braking cycle. These figures 
are contrasted to a total skid time of 4.47 sec. and a maximum magni- 
tude of skid of 68 ft /sec. for the uncorrected situation. 

When the brakes are initially applied at time = 0 sec. there is 
a time lapse of 0.16 sec. before a skid occurs. This is exactly the 
same situation that occurs in the uncorrected braking and is accounted 
for by the same reasoning. After 0.16 sec., a skid commences. This 
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Fig. 16. Double reap control response, K-j * 100000 lbs/sec., 
JC 4 = 100000 lbs/ sec. 



^■7 



is indicated in Fig. 15 by the positive going error signal and also 

• 2 

by the abrupt drop of the cu plot from - 24 rad/sec. to approximately 

2 

- 78 rad/sec. . Recall that when skidding occurs the net torque 

acting to slow the wheel is suddenly increased, thus the wheel 

starts to slow down rapidly, as seen from inspection of the cu 

versus time plot at the time of the initial occurrence of the first 

skid. However, the corrective force immediately begins acting to 

decrease the net force F (= F. - F ) acting on the brake shoe, 

which consequently decreases the net torque T (= ^ 2 ^ ~ 

and the wheel begins to speed up to "catch up" with the velocity of 

the vehicle. This accounts for the sudden positive slope of the cu 

plot almost immediately after an error signal occurs, and also for 

the upswing on the w versus time plot after the initial drop of cu 

with each occurrence of a skid. The corrective force increases until 

it equals F^ and/or e = 0 once more, at which time it begins decreasing 

towards zero until a skid occurs again. 

An expanded plot of (JO and v versus time for two skid cycles is 

shown in Fig. 17. To obtain this plot the equations of motion of the 

system were solved numerically by slide rule with the value of = 

, 191 

5000 lbs/sec. and = 8000 lbs/sec. Also, the simplifying assumption 

was made that the coefficient of friction between the brake drum and 

brake shoe was constant at 0.8. While these values for K 0 and K, 

5 4 

and the simplifying assumption for were not used in the analog 
computer runs, the shape of the curves is similar to the shape of 
those obtained on the analog computer simulation. A brief explanation 
of Fig. 17 will facilitate understanding the curves shown in Fig. 15 
and Fig. 16. 
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Fig. 17. Expanded plot of wheel Telocity and vehicular velocity 
versus tine for the double ranp control. 



Referring to Fig. 17, note that a skid begins, at time = 0 
sec. This is due to the simplifying assumption of a constant 
10. = 0.8 over the entire velocity range of interest, and in no way 
affects the qualitative aspects of the curve shown. The value of 
U) (or R^u) as plotted in Fig. 17) decreases initially and then rises 
to intersect the velocity V plot (which is decreasing linearly as 
the vehicle skids). For time less than the point marked (T) , the 
slope of the curve would be positive, as can be seen by mental 

• 

differentiation of the R^w curve. This is in agreement with the u) 

curve of Fig. 15 and Fig. 16. At the point marked (T) , however, 

the RjW) curve now begins increasing linearly, and the w curve would 

have a constant value, i.e., would level off at some positive value 

corresponding to the slope of the segment from (T) to Q) . This 

leveling off of the plot is observed also in the analog computer 

results of Fig. 15 and Fig. 16. The time period when the leveling 

off occurs corresponds to the region between (T) and (?) of the 

inset plot of F , i.e., where F has attained its maximum value 
r c c 

and the net torque is now a constant. Since the net torque is 
constant, the value of U) is also a constant, until the skid is 
completely eliminated. 

At point (2) the skid is eliminated and F^ begins dropping to 
permit reapplication of the full brake force F^. However, the net 
force F increases gradually from zero with a slope K^, and for a 
significant time after the elimination of the skid the velocity 
is decreasing at an unacceptably slow rate. This is seen by an 
inspection of the region of the curve of Fig. 17 from point (2) to 
(4) . At point (3) the deceleration has finally increased to that 
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value of deceleration experienced in a skid, i.e., - 16 ft/sec. , 
and at point ( 4 ) it is sufficiently great to cause a skid to occur 
once more. The process then repeats itself. 

The rather detailed qualitative analysis undertaken above serves 
to point up the major defects inherent in the double ramp form of 
control. If the parameters and are chosen so as to maximize 
the time when a skid is not occurring, e.g., to maximize the time 
from point ( 2 ) to point ( 2 ) of Fig* 17, then the time from point ( 2 ) 
to ( 3 ) will also be maximized. From ( 2 ) to (3) the velocity is 
decreasing at less than the optimum rate, and even at less than 
the deceleration experienced when in a skid. Thus it would seem 
that by minimizing the skid time with the double ramp control, we 
are inadvertently maximizing the actual time to stop. This conclusion 
is borne out by the results of Fig. 18 and Fig. 19, where the skid 
times and stop times for the various values of and simulated 
on the analog computer are plotted. 

Figure 20 is a plot of the magnitude of the skid observed versus 
K^> with fixed at 40000 Ib/sec. The magnitude of the skid after 
it has obtained a constant value is the value that is plotted here. 

It is seen that the magnitude of the skid is very dependent upon the 
value of K^, which is an expected result since the larger the value 
of the more rapidly is the skid corrected. This also accounts 

for the decrease in the skid time with increasing observed in Fig. 
18. Experimental results indicated that the magnitude of the skid is 
unaffected by the value of K^, which is natural since the downward 
going portion of the double ramp control is only acting when no skid 
is present. The stop time is increased, however, with increasing 
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or K^, since, as mentioned previously, the greater the proportion 
of the total stop time that is free from a skid, then the greater 
will be the actual time when the deceleration is less than that 
experienced when a skid is present. 

Even though the double ramp control has serious defects, it 
nevertheless represents a decided improvement over the completely 
uncontrolled braking represented in Fig. 14 0 Since under many road 
conditions it is the loss of control which occurs when a skid is 
present that is actually more dangerous than the increase in the 
stop time, then the double ramp control would have certain merits. 
The double ramp control is capable of minimizing the skid time and 
the magnitude of the skid, at a cost of slightly increased stopping 
distance „ 

If it is assumed that the stop time, the skid time, and the 
degree of skid experienced are of equal importance in judging the 
overall goodness of a particular type of control, then a performance 
factor (P.F.) can be defined for the uncorrected case, the double 
ramp control, and for subsequent controls to be tried. This per- 
formance factor is defined as the reciprocal of the product of the 
actual stop time minus the optimum stop time, the skid time, and 
the degree of skid experienced* The degree of skid will be measured 
after the velocity has decreased below 60 ft/sec., since that is the 
velocity below which the value for |JL (between the brake drum and 
brake shoe) is a constant. The optimum stop time is 3.13 sec. 

This definition of a performance factor must not be taken to 
be the absolute criterion for judging a control and can only be 
used as a rough guide to determine the relative merits of different 



45 



types of controls. It cannot, for example, be used to determine 
the proper setting of control parameters to give the ideal stopping 
conditions, since the optimum settings for one set of road conditions 
might very possibly be different from the optimum settings for another 
set of road conditions. Also, no account is taken, in the definition 
of a performance factor, of the fact that a slight degree of skidding 
might be much less significant than a larger degree of skidding, or 
that an increase in the stop time might be permissible under some 
circumstances but entirely unacceptable under others. Nevertheless, 
the definition given above is convenient and useful, within limits. 

For the uncorrected situation of Fig. 14, a performance factor 
can be computed in order to provide a basis for comparison. The 
integral of the error curve, rather than the product of the error 
times the skid time, is used, multiplied by the total stop time minus 
the optimum stop time. This gives a performance factor of 0.00214 
for the completely uncontrolled situation. 

For the double ramp control, the values of the performance 
factor range from 0.0154 to 0.298, and the variation of P.F. with 
a change in is shown in Fig. 21. This figure is presented only 
to indicate the trend of P.F. for the double ramp control and not as 
an aid to selecting the "best" value for . However, based on the 
indication given by the performance factors for the double ramp 
control, the double ramp control provides a definite overall perform- 
ance of the system, when compared to the uncontrolled situation. 
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6. THE BANG-BANG CONTROLS 0 



The bang-bang control with a dead zone was simulated on the 
analog computer. The dead zone was arbitrarily set at 7.5 ft/sec., 
which is rather high for a physical system, but which nevertheless 
serves to illustrate the behavior of the system under this form 
of control. The response obtained is shown in Fig. 22. 

The stop time is seen to be 4.50 sec., which is slightly less 
than that obtained for the completely uncontrolled situation. This 
is a predictable result, since for almost the entire cycle the 
vehicle is in a "controlled 11 skid, where the term "controlled" 
implies that the value of the skid is never permitted to exceed a 
certain value. The reason for the slightly lower stop time than 
that observed for the uncontrolled situation is that for a small, 
non-zero time interval in the region where F c drops to zero and F 
rises to F , the value of the skid is essentially zero (or at any 
rate less than the threshold established for detection of an error 
on the basic system of Fig. 7). While the skid is this small, the 
coefficient of friction controlling the deceleration of the vehicle 
is (which is greater than p^) , and the deceleration is thus 
greater for this small time interval. Since this phenomenon occurs 
every time the skid is corrected, the cumulative effect over the 
entire stopping distance is to decrease the stop time slightly. 

The major disadvantage to this form of control is the fact 
that a skid, even though of a small value, is almost always present 
until the vehicle comes to a complete stop. This also implies that • 
the shortest possible stop time is nearly 4.64 sec., i.e., the same 
time required to stop when no corrective force is present. However, 
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the fact that the skid is controlled, and can theoretically be 
reduced to an arbitrarily small positive value (by reducing the 
dead zone) indicates that even this form of control is superior to 
the uncontrolled situation. An additional advantage is that as 
the degree of skid is reduced toward 0, i.e., as the relative 
velocity between the tire and road surface is reduced, the coef- 
ficient of friction becomes somewhat greater than that present 
when a more significant relative velocity is attained. This means 
that a greater value of deceleration is possible when the skid 
becomes extremely small. This feature of the coefficient of 
friction was not incorporated in these analog computer simulations 
due to the greatly increased complexity of the simulation that would 
be necessitated. 

Even though this control could theoretically provide good perform- 
ance if the dead zone were reduced sufficiently, there is a serious 
physical limitation that prevents this from ever really being possible. 

A little reflection will reveal that the lower the dead zone setting, 
the more rapidly will the corrective force be required to go on and 
off. As the dead zone approaches zero, the shape of the corrective 
force will appear as a high speed oscillation, and the practical lower 
limit for the dead zone will be determined by both the speed of response 
of the hydraulic or pneumatic lines and the physical strength of the 
brake system. 

Thus, from a practical point of view, this system appears to offer 
the advantage of a controlled skid, plus a shorter stop time than that 
provided by the double ramp control. The disadvantage is that a skid 
is almost always present, and the level to which it could be reduced 
would be severely limited by the physical response of the brake lines. 
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Neither control tried thus far appears to be superior to the 
other* This is verified somewhat by determining the performance 
factor for the bang-bang control with a dead zone. The performance 
factor (computed by taking the reciprocal of the product of the 
integral of the error and the difference between the stop time and 
the optimum time to stop) is determined to be 0.046. This is 
neither significantly better nor worse than any of the performance 
factors obtained for the double ramp control. 

The bang-bang control with a time delay was then simulated. 

The results for this simulation are presented in Fig. 23. A time 
delay of 0.1 sec. (problem time) was chosen for this simulation. 

The stop time is seen to be 5.93 sec., which is considerably worse 
than the stop time for the previously described bang-bang control. 
However, in this case the skid is eliminated for a portion of the 
total stop cycle, and the skid time is 3.92 sec „ Superficially this 
appears to be an improvement over the previous bang-bang control, 
but closer examination reveals that very little, if any, betterment 
of the performance is gained by the incorporation of a time delay. 

Reference to Fig. 24, which is an expansion of a portion of 
Fig. 23, and was obtained by manual solution of simplified equations 
of motion, reveals that for a significant period of time during each 
error correcting cycle the velocity remains a constant. For example, 
from point (Y) to (J) , where no skid is present, the velocity is 

also a constant. The time from (j2) to Q) is 0.1 sec., and corresponds 
to the time delay built into the corrective force. During this time, 
the net force F on the brake shoe is zero, and thus the velocity 
remains constant. If the time delay is decreased, this control 
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Fig. 24. Expanded plot of wheel Telocity end Tehicular Telocity 



approaches the performance of the bang-bang control with a zero 
dead zone. 

The bang-bang control with a time delay thus appears superior 
to a completely uncorrected braking situation, but appears some- 
what inferior to either of the two controls previously described. 

A performance factor for this control (for the specific set of 
parameters used for this simulation) is 0.025, which is sufficiently 
lower than the trend of the values obtained for the previous two 
controls to warrant the above conclusion. 
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7. THE TRIPLE RAMP CONTROL . 



The triple ramp control of Fig. 12 was the final type of 
corrective scheme to be implemented a A total of 100 separate 
computer runs were made, using a variety of values for the 
parameters K^, K^, and C.R. (correction ratio). Representative 
samples of the behavior to be expected from this type of corrective 
scheme are shown in Fig. 25 and Fig. 26. 

The most striking features about both Fig. 25 and Fig. 26 
are the low stop time, skid time, and magnitude of error present 
for this form of control. For example, Fig. 25, with = 60000 
lbs/sec., = 40000 lbs/sec., = 1000 lbs/sec. and C.R. = 2.05 
has a stop time of 4.22 sec., skid time of 0.98 sec. and roughly 
1.9 ft/sec. as the degree of skid. This gives a performance factor 
of 0.493 which is higher than the maximum value attained for any 
control tried previously. This value of performance factor is of 
course affected by the values of the various parameters K^, K^, , 

and C.R., and does not represent the "best" value of P.F. that was 
obtained on the trial runs for the triple ramp control. Indeed, 
the term "best" is misleading here, since a setting of parameters 
resulting in the highest possible P.F. would sacrifice an adapta- 
bility to changing road conditions. This feature will be discussed 
more thoroughly later in this paper. 

Nevertheless, it can be concluded from the above factors and 
an inspection of Fig. 25 and Fig. 26 (along with the Author* s 
analysis of the other 98 analog computer runs) that the triple ramp 
control provides more nearly optimal braking conditions than any 
previously described control. 
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A closer examination of the actual behavior of this triple 

ramp control is appropriate at this point. The results of a 

manual solution of the equations of motion (with the simplifying 

assumptions mentioned previously for similar hand calculations) are 

presented graphically in Fig. 27. The parameter values assumed 

were K 0 = 50000 lbs/sec., K, = 80000 lbs/sec., K = 1000 lbs/sec., 

3 4 c 

and C.R. = 1.5. Also, the coefficient of friction |JL between the 
brake drum and brake shoe was assumed to be constant at 0.8 for 
the hand calculation (this was the only simplifying assumption). 

All other conditions were assumed the same as previously listed 
for other simulations. 

An examination of Fig. 27 indicates that from time = 0 sec. 
to the point marked (3) , the behavior of the control and the 

vehicle response are precisely analogous to the double ramp form 
of control. This can best be seen by a comparison of Fig. 27 
with Fig. 17 over the region from time = 0 sec. to the point ( 3 ) . 

The same defect present in the double ramp control is also present 
in the triple ramp control, i.e., from point ( 3 ) to point ( 5 ) 
the velocity is decreasing at a rate slower than the rate of 
deceleration when a skid is present. Also, from point ( 3 ) to point 
( 6 ) this phenomenon occurs. However, at point (j) , when the down- 

ward slope of the corrective force F c abruptly changes from - 
to - the situation changes completely. While the slope of 
F^ = - K^, the rate of decrease of the velocity not only is signif- 
icantly greater than that possible during a skid, but also no error 
is present until the value of the deceleration again exceeds - 24.0 

2 I—. 

ft/sec. at point (j5 ) . At this time the process repeats itself. 
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▼ersus time for the triple reap control. 



It is during this last phase of the skid correcting cycle 
(from point ( 7 ) to point (8) in Fig. 27) that the inherent 
superiority of the triple ramp control is manifested. During 
this portion of the cycle the deceleration is gradually increasing 
as F^ slowly decreases. When the deceleration is too great for 
the prevailing road conditions a skid once again occurs, but prior 
to this occurrence of the skid the deceleration is more and more 
closely approaching the optimum deceleration. Additionally, if 
the road conditions are changing, for example, if the vehicle 
went from a patch of slick pavement to dry pavement, and the coef- 
ficient of friction between the tire and road surface were increased, 
then this form of control permits the deceleration to increase to 
the new optimum value. In other words, the triple ramp control 
is continually sampling the road conditions and permitting the 
braking force to adapt to changing road conditions. 

The desirability of this adaptability feature is evident. It 
is also apparent that for the invariant road conditions simulated 
for this problem, the most nearly ideal braking situation would 
have been achieved by setting C.R. to a value just slightly greater 
than one, and causing to equal zero. This would have permitted 
only two skids of small magnitude to occur, and for the remainder 
of the path the vehicle would have slowed down at a rate very nearly 
equal to the optimum, deceleration. This would provide a very large 
value for the performance factor, which would of course be a misleading 
indication in this case, since the control would not respond to improved 
roaxl conditions. It would, however, still respond to worsening road 
conditions, so there is even some merit in this proposal. 
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Referring back to Fig. 25, we note that the typical triple 
ramp shape of F^ does not become evident until after the occurrence 
of the fourth error signal, whereas in Fig. 26 the triple ramp 
shape appears as expected after the occurrence of the second error 
signal. Recall that the form of the coefficient of friction p. 
between the brake drum and brake shoe is given by 

= .8 - .01 (u) - 60) for uo > 60 rad/sec. 

Reference to the uo versus time plot of Fig. 25 reveals that both 

the first and second errors occur while the rotational velocity is 

above 60 rad/sec. Thus in this region the value for p* is constantly 

changing in accordance with the above relation. Now, when the first 

error signal (or skid) is sensed the corrective force F^ immediately 

increases until the skid is eliminated, and then decreases until 

the skid occurs once more: this value of F at which the skid occurs 

5 c 

is stored by a storage device. Upon the occurrence of the second 
skid, F^ increases again until the skid is once more eliminated. 

Ideally the value of F^ should drop to some multiple of the stored 
previous value of F^ (the multiple being established by C.R.). How- 
ever, the velocity of the vehicle is decreasing continually during 
this process and p. is thus increasing. This means that after the 
occurrence of the second skid, the net torque on the wheel for a 
given value of F^ (e.g., for a value equal to C.R. times the previously 
stored value of F ) is now greater than it would be if the value for p. 
were not increasing. Thus a third skid occurs before F c can drop to 
the value at which the break point occurs. However, after the occur- 
rence of the fourth skid, p. is constant, and thus the typical triple 
ramp effect is observed. 
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In Fig. 26 the above phenomenon does not occur. Note that 

in this figure the value of is only 25000 lbs/sec. whereas in 

Fig. 25 had a value of 60000 lbs/sec. This means that the 

error will be corrected more slowly, and thus the velocity will 

have decreased to a level where is a constant by the time the 

second error has been eliminated. Thus the coefficient of friction 

will not have changed sufficiently between the time when the value 

of F was stored and the time when F reaches its break point to 
c c r 

cause another error signal to occur. Note, however, that if C.R. 
had been set sufficiently low (e.g., slightly above 1.0) then the 
same phenomenon that occurred in Fig. 25 would also have occurred 
in Fig. 26. Conversely, had C.R. been set sufficiently high, this 
phenomenon would not have occurred in Fig. 25. 

Figure 28 presents a plot of stop time versus , with C.R. 
as a parameter. The data points on this and subsequent plots 
were obtained from the graphical output resulting from the various 
analog computer runs made while simulating the triple ramp control. 
Consequently, the accuracy of the individual data points is only 
of the order of a few percent. However, this degree of accuracy is 
sufficient to demonstrate the overall performance characteristics 
of the triple ramp control quite satisfactorily. From Fig. 28 we 
can see that as K increases the stop time increases also. However, 
the rate of increase, i.e., the slope of the curve, is observed to 
be a variable, even though the overall trend is upward. This over- 
all upward trend is expected in view of previous comments concerning 
the adaptability feature inherent in the triple ramp control. A 
decrease in stop time must be achieved at the cost of a decrease in 
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the sensitivity of the control to improving road conditions, all 

other parameters being assumed constant. Thus, as increases, 

the vehicle goes into a skid mode more frequently, and the overall 

trend of stop time with K would be an upward one. 
r c 

However, to satisfactorily account for the leveling off of 
the stop time versus chart at localized regions, we must investi- 
gate the detailed behavior more closely. Figure 29 is a plot of 
both stop time and skid time versus for one specific set of para- 
meters. The numbers adjacent to the skid time curve represent the 
total number of skids present in the particular graphical output 
from which each specific data point was obtained. For example, 8 
skids were present in the analog computer output obtained with a 
setting of = 60000 lbs/sec., = 25000 lbs/sec., C.R. = 2.05, 

and K = 1500 lbs/sec. For K = 2500 lbs/sec. and K = 3000 lbs/sec., 
c c c 

10 skids were present and the skid time was thus the same for each 

of these values. However, Fig. 29 also reveals that the stop times 

for each of these values of K were the same. 

c 

In a situation where the same number of skids is observed from 

one value of K to another, then the skid times will obviously be 
c 

identical, assuming that (which directly determines the magnitude 
of each skid and the length of time it is present) is the same for 
each trial. The only factor contributing to a difference in over- 
all stop time is thus the difference in time between adjacent pairs 
of skids (e.g., the time from point (i) to point (8) in Fig. ^ 27). 

For a larger value of , the time between pairs of skids is decreased 
Now, for example, when 10 skids occur for a value of K 3000 lbs/sec 
the time between adjacent pairs of skids is shorter than the time 
between the adjacent pairs of 10 total skids obtained for = 2500 
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lbs/sec. Thus the final pair of skids occurs earlier when K = 

c 

3000 lbs/sec., but due to a larger value of a larger value of 

deceleration is attained after the tenth skid and prior to the 

vehicle coming to a complete stop. There is also a longer period 

of time after the tenth skid occurs, during which the slope of 

F = - K , for the case where K = 3000 Ibs/sec. These two 
c c c 

factors serve to account for the leveling off of the stop time 
versus curves observed for Fig. 28 and Fig. 29. It is also 
possible that a slight local decrease in total stop time might 
occur for an increase in under the proper combination of 
control parameters and road conditions. This is observed in Fig. 

28 in one region of the curve obtained for C.R. = 2.74. 

Despite the above observations, the conclusion to be drawn 
from Fig. 28 and Fig. 29 is that the overall effect of an increase 
in is to increase the stop time of the vehicle, and to degrade 
the overall performance of the corrective system. A compromise 
must thus be reached between adaptability to improving road 
conditions and decreasing the total stop time. Designing for a 
worst case situation, a satisfactory value for would thus be 
= 500 lbs/sec., as this would permit reasonably rapid response 
to improving road conditions while still providing very good braking 
performance under nearly static road conditions. 

The next parameter to be investigated is the correction ratio 
(C.R.). Figure 30 is a plot of the stop time versus the correction 
ratio, with as a parameter. The overall trend as determined from 
Fig. 30 is that the stop time tends to decrease with an increase in 
the value of C.R. Each curve in the family of Fig. 30 demonstrates 
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notable local exceptions to this trend, however, and, as will 
be pointed out, excessively large values of C.R. will materially 
lengthen the stop time. 

The overall downward trend of the stop time versus C.R. curve 

is best understood by referring to the inset of Fig. 27, where a 

pictorial description of the triple ramp control is presented. 

As the value of C.R. is raised this raises the break point at 

which the transition of the slope of from - to - occurs. 

This has the effect of lengthening the time between the occurrence 

of adjacent skid pairs, since the time that F^ will take to drop 

to the level where a skid will occur once more increases in direct 

proportion to the height of the break point. It is thus seen that 

increasing C.R. has the same net effect as decreasing the value 

of K c , with certain minor restrictions. 

It would be expected that, as C.R. is increased to obtain a 

set of braking histories, those output traces which have the same 

number of skids present for differing values of C.R. would behave 

in a manner similar to that observed when K was the variable 

c 

quantity (e.g., in Fig. 28 and Fig. 29). Referring in Fig. 30 to 

the curve obtained for K = 500 lbs/sec., we see that the stop time 

c. 

increases with increasing C.R. when the total number of skids present 
in each braking history remains a constant four (this is also observed 
for = 1000 lbs /sec., for six skids observed). The reason for 
this is very simple. As C.R. is raised this "pushes" the adjacent 
pairs of skids farther apart and thus for a larger portion of the 
total stop time the vehicle is decelerating at a lower rate than if 
C.R. were set lower. Thus if the total of skids still remained the 
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same, the control which provided the greater deceleration over 
the entire skid free path would also provide a shorter stop time. 
This is the control with the lower C.R. 

Note that C.R. is also closely linked to the adaptability of 
the system to improving road conditions, since for a large value 
of C.R. the time required for F^ to drop (at a rate of - K^) to 
a level where the road is ’’tested 11 for better conditions is 
lengthened. Additionally, C.R. cannot be increased to a very 
large value without greatly extending the stop time. The minimum 
number of skids possible (starting with any significant initial 
velocity) is two. Once C.R. has been increased to a level where 
only two skids are present, then any further increases in C.R. will 
cause a rapid increase in the stop time for the reasons listed in 
the preceding paragraph. However, long before the skids were 
reduced to two in number, the stop time would also have increased 
due to the low value of deceleration that would be present over an 
excessively large portion of the vehicle’s trajectory. 

In retrospect, it appears that despite the downward trend of 
stop time with increasing C.R. presented in Fig. 30, any advantages 
to be gained by very large values of C.R. are outweighed by the 
disadvantages. On the basis of information presented here, a ’’safe” 
value for C.R. would appear to be about 2.0. 

Figure 31 and Fig. 32 present plots of both stop time and skid 
time versus C.R. for two different sets of parameters. They are 
included primarily for information, the most significant factor 
about either one being that where the curve representing skid time 
becomes flat (implying a constant number of skids were observed for 
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Fig. 32. Stop tine and skid tin* Tarsus C.R. for the 
triple rasp control, * 60000 lbs/ sec., 

* 40000 lbs/sec., K c * 1000 lbs/ sec. 
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succeeding braking history plots) the curve for stop time acquires 
a positive slope. This serves to reinforce the arguments previously 
advanced. 

The behavior of the braking performance as is varied is 
also of some interest. Figure 33 illustrates the behavior of stop 
time and skid time as is increased. Notice that the tendency 
exhibited is for the stop time to decrease with increasing . This 
is the expected result, since the larger the value of , the more 
rapidly is any error or skid corrected back to zero, and thus a 
larger portion of the braking path will be in a controlled mode of 
skid free braking. 

Referring back to Fig. 27 once more, we note that the time 
from the commencement of braking action (t = 0 + ) to point (T) , and 
also from point (4) to point (3T) , is dependent solely upon the value 

of K^, whereas the time from point (2) to point (4) and from point ( 5 ) 
to point (8) depends upon and . Thus as is increased the 

total skid time and magnitude of each individual skid are reduced. 
This has the effect of causing each successive skid pair to occur 
earlier in time, even though the time between skids is unchanged. 

As long as the total number of skids present is unchanged as is 
increased, an improvement will always be observed in the stop time. 
However, in the event that one value of causes a certain number 
of skids, and a higher value of is just sufficient to increase 
the number of skids, then a slight local increase will be observed 
in the stop time. This phenomenon is seen to occur in Fig. 33. 

The explanation for this is as follows. When = 50000 lbs/sec. 
the time remaining between the correction of the sixth skid and the 
complete stopping of the vehicle is insufficient for the deceleration 
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SKID TIME (Seconds} 



to have increased to that level necessary to cause another skid. 

However, for = 60000 Ibs/sec. the seventh skid occurs just 

prior to the vehicle coming to a rest. During this seventh skid 

2 

the vehicle deceleration is only - 16 ft/sec. , rather than the 

2 

- 24 ft/sec. that was experienced just prior to the skid. It is 
this lower value of deceleration present just prior to stopping 
that causes the stop time to increase as is increased from 
50000 lbs/sec. to 60000 lbs/sec. 

The above phenomenon is just restricted to the transition 
points where the number of complete skids increases from one value 
to another. However, the decreased skid times and magnitude of 
skids caused by an increased value of indicates that the greatest 
overall advantage is to be gained by making as large as possible, 
regardless of any localized increases in stop time caused by this 
policy. 

Finally, the system behavior with changing values of is 
investigated. Figure 34 indicates that the stop time decreases as 
the value of is increased, but that the skid time goes up with 
increasing K^. Reference to Fig. 27 again facilitates understanding 
this behavior. From Fig. 27 it can be seen that the time lapse from 
point (T) to point (4) and also from point (3) to point (j) is depend- 
ent solely upon the value of K^. Large values of will cause the 
time intervals between the points to become small, while leaving 
time intervals along other portions of the trajectory unchanged. 

Since it is during the time interval from point (5) to point (§) , 
for example, that the deceleration of the vehicle is even less than 
that present while in a skid, any factor which will shorten the time 
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from (5) to ( 7 ) will also shorten the time the vehicle spends in 
the undesirable region from (5) to (6) . This will naturally 
cause a shorter stop time, since the vehicle will be decelerating 
at a high rate over a larger proportion of its stopping path. 

The increase of skid time noted in Fig. 34 is once again 
related closely to the total number of skids present. As the time 
from (T) to (4) or from (5) to (7) (referring once again to Fig. 27) 
is reduced, this tends to shorten the time between the occurrence 
of adjacent skids. If the total number of skids remains the same 
while is being increased, the skid time naturally remains the 
same (to fully appreciate this, one must recall that the duration 
and magnitude of each individual skid is controlled by only, 
which is not being varied here). However, as the skid occurrences 
become closer together in time, a point is reached where the total 
number of skids encountered in the stopping history increases, and 
thus the skid time will also increase. This increase in skid time 
would not be a serious consequence, however, if were sufficiently 
high to make the duration and magnitude of the individual skids very 
sma 1 1 . 

Thus it might be concluded that the larger the value of K^, 
the better the performance of the entire system, since the decrease 
in stop time would more than offset any increase in skid time, partic 
ularly if were also set very high. Caution must be exercised here 
however, because it is the nature of the triple ramp control that 
the value of at which every second error occurs must be marked 
and stored. This means that if is too large, then difficulties 
would be encountered in marking the point at which each second error 
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occurs, and additionally in reversing the trend of (i.e., in 
changing from a slope of - to + Obviously a slope of 

- infinity for - would be totally unsatisfactory. Thus, a 
very large value of would indeed provide the best performance, 
but this value must not be so great as to create excessive dif- 
ficulties in measuring and storing the point at which each second 
error occurs. 

In summary, the following qualitative comments concerning 
the choice of the values for the various parameters of the triple 
ramp control are appropriate. should, with no qualifications, 

be made as large as possible. should also be as large as possible, 

short of being a negative step function, and subject to the restric- 
tions mentioned in the preceding paragraph. K c should be relatively 
low assuming that C.R. is not too large. A value for C.R. of about 
2.0, and for of 500 lbs/sec. was determined to be satisfactory 
for the particular system simulated here. 
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8. PROPOSED IMPLEMENTATION OF THE TRIPLE RAMP CONTROL „ 



In view of the obvious superiority of the triple ramp form 

of control over any of the other controls examined, a proposal 

for the physical implementation of this control is advanced here. 

It must be assumed that a linear accelerometer of high accuracy 

2 2 

and linearity over the range from 0 ft/sec. to 32 ft/sec. is 
available,, A high gain, broad bandwidth operational amplifier is 
necessary to provide the integrator used in conjunction with the 
linear accelerometer (see Fig. 1 for the block diagram of the 
system). Additionally, an accurate tachometer to measure the 
rotational velocity of the wheels is assumed to be available. 

Thus, all elements necessary to accurately register the occurrence 
of a skid are assumed to be available. Our attention will be 
directed toward the actual mechanism by which the triple ramp 
control is generated, and by which the pressure applied to the 
brake cylinder is modified. 

Figure 35 illustrates the means by which it is proposed to 
implement the triple ramp control. An explanation of Fig. 35 
follows . 

A force which is applied to the brake pedal is transmitted to 

the brake drum via the hydraulic lines. Until a skidding situation 

occurs is held equal to V ^ by means of the ramp generator and 

limiter. When V equals V , full braking action is applied and 
c max 

no corrective action is occurring. 

The first time a skid occurs, e (error) becomes greater than 
zero, and switches (T) and (T) throw in the indicated manner. 
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Capacitor discharges rapidly through resistor (a very small 
resistance). thus drops very rapidly, permitting the braking 

force to be decreased toward zero; this step is an exponential 
decay with a very small time constant, but it approximates the 
first upward ramp of (shown with a slope of K^). 

When the error is reduced to zero (skid eliminated) switches 
(l) and (2) throw once more. The voltage on then starts to climb 
with a slope of K^, thus causing the braking force to increase once 

more. During this first part of the cycle a portion of V c is being 

monitored by (thru the tap on R^) . When the braking force has 
increased sufficiently to cause skidding to occur once again, 
switches (l) and (2) throw once more and drops rapidly, decreasing 
the brake force as before. Additionally, however, switches (J) and 
© are now thrown, and capacitor has a charge proportional to 
the force required to cause skidding (say . 9 x to cause a skid). 

When falls sufficiently to eliminate the skid, e = 0 again, 

and the ramp generator charges upward once again. When is 

just equal to the charge on C- (.9 x V to cause a skid, for example) 

1. c 

V t t> ecomes greater than zero, and switch (jT) throws as indicated. 
This puts a large capacitance (C^) in parallel with and causes 
the ramp slope to change to a much shallower slope (i.e., to change 
from slope = to slope = K^) . is thus caused to approach the 

value that creates a skid much more slowly. 

When (or if) a skid occurs once more, the entire process 
described above repeats itself. 

A sketch of the force acting on the brake cylinder is shown in 
Fig. 36. 
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Fig. 36. Force acting on the brake cylinder during occurrence 
of one skid pair, if using Proposed Implementation. 

Note also that the sensing of the first, third, etc., skids 
or second, fourth, etc., skids can be accomplished by an R - S 
flipflop connected to the output of a comparator which would be 
used to indicate the presence or absence of a skid. This was 
done on the analog computer simulation with very satisfactory 
results . 

Figure 37 is a blow up of the pressure controller shown in 
the inset of Fig. 35. This is a standard form of control valve, 
and is not original with the Author. [l] A brief explanation of 
the workings of a standard power brake system, for which this 
particular form of implementation is ideally suited, is appropriate 
here . 

In a power brake unit, when the brakes are applied, a valve 
to the atmosphere is opened on one side of a piston, and the other 
side is connected to a vacuum. The atmospheric pressure forces 
the piston closed, and causes pressure to be transmitted to the 
brake cylinder via a hydraulic line. 
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Fig. 37. Detail of pressure controller suitable for use 
with a power brake system. 



The above controller will cause a disruption of the pressures 
on the power cylinder if the valve is rotated from the position 
shown. In the position shown, the vacuum line is blocked off thus 
permitting full atmospheric pressure to the piston; any other 
position of the valve would lessen the pressure on the piston, 
thus reducing the braking action. 
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9. CONCLUSIONS. 



Of the three major types of controls simulated, the triple 
ramp control proved to be the only one lacking major defects, in 
that it reduced the degree of skid to a very low level while still 
providing a very short stop time. Nevertheless, the double ramp 
control was determined to provide a definite improvement over the 
completely uncontrolled situation, as did the bang-bang controls. 
However, for the double ramp control a slightly increased stopping 
time was the penalty paid for obtaining a lower degree of skid 
while stopping the vehicle. The bang-bang control with a dead 
zone proved unsuitable due to the presence of a slight degree of 
skid during the entire braking path of the vehicle. Additionally, 
the level to which the skid could be reduced would be severely 
limited by the response characteristics of the brake lines. The 
bang-bang control with a time delay, while still superior to the 
uncorrected situation, appeared to be considerably inferior to any 
other form of control simulated. 

The triple ramp control, however, appears to provide very 
nearly optimum braking conditions in that the degree of skid and 
the skid time are both minimized, while still providing a considerably 
shorter stop time than would be obtained from any other form of 
control. A trade off between stop time and adaptability to improving 
road conditions is necessary though, and the proper balance between 
these two factors would be determined by both the type of vehicle 
involved and the anticipated road condition. However, greater 
emphasis should probably be placed on minimizing skid time and stop 
time by having a relatively low value for K^. 
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The problem of possible oscillation introduced by the triple 
ramp control should never arise, since it is essentially a negative 

C6] 

feedback system and should be inherently very stable. Extremely 
unusual road conditions would probably be necessary for any oscil- 
lations to occur; however, an investigation of this subject has 
not been undertaken. This would possibly be a fruitful area for 
further study. 

The response time required for the mechanical components of 
the proposed implementation of the triple ramp control (Fig. 35 and 
Fig. 37) would certainly be easily attainable. This may be seen 
from Fig. 25 (which is a "typical" response curve), where examina- 
tion of the F c versus time curve indicates that the pressure control- 
ler would go through the high speed cycle at the rate of only about 
4 c.p.s. This serves to give an order of magnitude approximation 
to the response characteristic necessary. 

Throughout this entire thesis, F has been assumed to be a 

constant, implying a panic stop situation with maximum force applied 

to the brake pedal. If, during the course of slowing the vehicle, 

the operator decreased the pressure applied to the brake pedal, no 

effect would be noticed in the operation of the automatic control, 

as long as the value of F were still sufficiently large to cause 

a skid (if acting by itself). If, however, F^ n were decreased below 

the value needed to cause a skid (if acting alone), then the control 

would permit the net force F to rise to the new value of F. (either 

in 

at the rate of or K^, depending upon which part of the cycle was 
acting at the instant of decreasing the value of F^) and this value 
of F would be retained until a skid once again occurred. If the 
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applied brake force were increased, the action of the triple 

ramp control would be unaffected, as the control would still 
continue to seek the optimum net force F. 

A matter of possibly more concern is the characteristic of 
the linear accelerometer and of the integrator used to indicate 
the existence of a skid. The proper operation of this skid detect- 
ing scheme is basic to the operation of any skid correcting device 
proposed in this thesis. The mounting of the linear accelerometer 
would be very critical, and a position on the axle housing, in 
order to minimize "nose dive" effects, would be the most nearly 
ideal location. However, no serious theoretical limitations on 
the effectiveness of this form of skid detecting scheme present 
themselves, and this would appear to be the only practical means 
of detecting a skid prior to actual wheel lock up. 

A final comment concerning the analog computer on which the 
brake system and all controls were simulated is appropriate at 
this time. The COMCOR Ci 5000 Analog Computer is in actuality a 
form of hybrid computer, in that extensive logic circuitry is avail- 
able for use in conjunction with standard integrators and operational 
amplifiers. Simulation of this type of problem would be difficult 
if not impossible on a more conventional analog computer lacking 
the logic circuitry, high speed switches and overall versatility 
of this particular computer. A simulation of the type described 
would be possible on a digital computer utilizing one of the 
languages available for this purpose, such as DSL-360, etc. How- 
ever, in order to obtain the extensive data needed to determine the 
performance characteristics of the various controls, large quantities 
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of very expensive digital computer time would be required. On 
the analog computer, however, almost two hundred different runs 
were made consecutively in a relatively short (and inexpensive) 
time period by merely varying potentiometer settings between runs. 
In this case, the accuracy of a digital computer was unnecessary, 
and the (hybrid) analog computer proved to be an exceedingly 
valuable design aid. This thesis thus illustrates graphically 
the continuing need for increasingly more sophisticated analog 
computers, such as the COMCOR Ci 5000. 
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APPENDIX A 



Determination of the maximum possible deceleration of a vehicle. 

The maximum deceleration possible for a vehicle is determined 
by the coefficient of friction between the tires and the road 
surface. Figure A1 represents an idealized vehicle traveling 

o • • 

in the plus direction with a velocity X, and a deceleration X, as 
indicated. 



»• 




Fig. Al. Idealized vehicle. 

In Figure Al, F^ represents the force of friction acting on 
the vehicle, and is the only external force causing the vehicle 
to decelerate. The weight, mg, is the normal force acting on 
the ground, so the value of F^ is given by 

F f = p*ng ) (Al) 

where is the coefficient of friction between the wheel and the 
road surface. The coefficient of friction is a variable, which 
has its greatest value (p, g ) just before a skid occurs, and then 
drops abruptly to a lower value p* when a skid begins. 

IX 

Thus the maximum value of F^ is 

F f (max) " 'V 8 • (A2) 
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Note that ^ corresponds to the coefficient of kinetic friction, 
whereas corresponds to static friction. 

Summing the forces in the X direction, we obtain 



mX = - F 



f * 



(A3) 



Substituting Eq. (A2) in Eq . (A3) and solving for > we 

obtain for the maximum possible deceleration 



" X (max) " 



(A4) 
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APPENDIX B 



Exp l a n ation of —log computer symbols used. 

Tho analog computer elements used in the simulations described 
in this thesis are indicated below. 



SYMBOL 




X, 




Xou-t 







X-re«T > K 

— y— 



X 



«u-f 




DESCRIPTION 

Sign o hanging summer amplifier. 

Xout - ~( x i +■ I OXj) 



Sign changing integrator. 



r 



Xout 



J' (x,+-iox £ ) dt +x c 



Class HI multiplier. 

Xout = x,x 2 

Potentiometer . 

Xout «= K| X| , K, < I 

Double pole, double throw relay. 

Xout = X| ) "for X Te sT < K 

Xout “ ) "for Xtest K 

High speed digital/ analog switoh. 

Xout = 0 ;-fot X Ttf « T <K 

Xout — X \ ) -f-of- XrecT> K 

Comparator. 

C » logical 1 and C • logical 0, 
whom xi + *2>0. 

C * logical 0 and C * logical 1, 
when xj + X 2 < 0 . 

Comparators are used to throw both 
the DPDT relays and the D/A switches. 
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TFST>>Kj 






Ti® point. 

Used in conjunction with two or wore 
D/A switch® a which sro thrown so that 

only on® switch is ON at a tin®. 

*0«t = x, or 

Switching acoonplish®d by MAMD logic, 
utilising only on® D/A switch. 

Xoat = X, j when 
CT5ST,>K,) and (TEST £ >K 2 ). 

Also , X 0 ut. = 0 , when 

(TEST, <K,) or ("tE STa< Ka) . 
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